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ABSTRACT 

We propose a "two-stage" model for the effects of feedback from a bright quasar on the cold 
gas in a galaxy. It is difficult for winds or other forms of feedback from near the accretion disk 
to directly impact (let alone blow out of the galaxy) dense molecular clouds at ~kpc. But if 
such feedback can drive a weak wind or outflow in the hot, diffuse ISM (a relatively "easy" 
task), then in the wake of such an outflow passing over a cold cloud, a combination of in- 
stabilities and simple pressure gradients will drive the cloud material to effectively expand in 
the direction perpendicular to the incident outflow. This shredding/expansion (and the corre- 
sponding decrease in density) may alone be enough to substantially suppress star formation in 
the host. Moreover, such expansion, by even a relatively small factor, dramatically increases 
the effective cross section of the cloud material and makes it much more susceptible to both 
ionization and momentum coupling from absorption of the incident quasar radiation field. We 
show that even a moderate effect of this nature can dramatically alter the ability of clouds at 
large radii to be fully ionized and driven into a secondary outflow by radiation pressure. Since 
the amount of momentum and volume which can be ionized by observed quasar radiation 
field is more than sufficient to affect the entire cold gas supply once it has been altered in 
this manner (and the "initial" feedback need only initiate a moderate wind in the low-density 
hot gas), this reduces by an order of magnitude the required energy budget for feedback to 
affect a host galaxy. Instead of ~ 5% of the radiated energy 100% momentum) needed if 
the initial feedback must directly heat or "blow out" the galactic gas, if only 0.5% of the 
luminosity 10% momentum) can couple to drive the initial hot outflow, this mechanism 
could be efficient. This amounts to hot gas outflow rates from near the accretion disk of only 
^ 5 — 10% of the BH accretion rate. 

Key words: quasars: general — galaxies: active — galaxies: evolution — cosmology: theory 



1 INTRODUCTION 

Observations have established that the masses of supermassive 
black holes (BHs) are tightly correlated with various host galaxy 
properties (Magorrian et al. 'l998'; 'perrarese & Merritt' 2000'; 
Gebhardt et al. 2000; Hopkins et al. 2007a; Alter c& Richstone 
2007.) . Together with constraints indicating that most of the BH 



feedback begins to expel nearby gas and dust. This "blowout" re- 
sults in a short-lived, bright optical quasar that, having expelled 
its fuel supply, fades and leaves a remnant on the observed BH- 
host correlations (Hopkins et al. 2005a, c). These scenarios have 



mass is assembled i n optically bright quas a rl phases i Soltanll 19821 ; 
ISalucci et alj|l999l ; IYu & Tremaind I2OO2I ; Hopkins et al.ll2006bl) . 
this has led to the development of models where feedback pro- 
cesses from accret i on self-regulated BH growth at a critical mass 



nosity functions, lifetimes, and BH mass functions ( 


Hopkins et al. 


2005b 


, I2OO6C. 2008bl, 2009 


; Volonteri et al. 


2006 


; Menci et al. 


2OO3I; 


Somerville et al.ll2008l; 


Lapi et al.ll2006l; 


Tortora et al.ll2009h. 



dSilk & Reesll998l ; lDi Matteo et al.l200^ ; lMurrav et alj2005l) . Gas 

inflows triggered by some process fuel rapid BH growth, until 
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' In this paper, we use the term "quasar" loosely as a proxy for high- 
Eddington ratio activity, rather than as a reference to specific optical prop- 
erties. 



It is much less clear, however, what the impact of whatever 
feedback regulates BH growth will be on the host galaxy. In mod- 
els, such feedback is invoked to explain the rapid "quenching" of 
star formation and sustained lack of cooling in massive galaxies 
(Granato et al. 2004; Scannapieco & Oh 2004; Croton et al. 200^; 
[Hopkins et al.i2008al ; lAntonuccio-Delogu & Silk 2008). The argu- 
ment in the models is that, under various simple assumptions, if 
sufficient energy or momentum is injected into the ISM near the 
BH on a timescale short enough to halt accretion, then it will yield 
a supersonic pressure or momentum-driven outflow that propagates 
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to large scales (see e.g . iMonaco & Fontano^l2005l : iHopkins et al.l 
l2006a| - |Shin e"tai]|2009t) . 

But the actual mechanisms of feedback and physics of the ISM 
rlevant for this remain highly uncertain. Highly ene rgetic outflows 
are as sociated with bright quasars (for a review, see lVeilleux et al.l 
l2005h : these range from intense winds (v ~ 10'*kms ') asso- 
ciated with the central engine seen in the broa d emission line 
regio ns and broad absorption line quasars (e.g. IWevmann et al.l 
Il98lh to more moderate outflows (v ~ 10^ — 10^^ km s~') associated 
with the narrow line region and the "warm absorber" jLaor et al] 
Il997l ; lcrenshaw et al.ll200ol) as well a s with small-scale quasar ab- 
sorption and occu l tation systems (e.g . iMcKernan & Yagoobll 19981 ; 
iTumer et alJl2008l ; lMiller et alfeOOsh . Indeed, high- velocity winds 
driv en near the accretion disk are theoretically hard to avoid (see 
e.g. iBlandford & Payne 'l982'; 'Beg elmaj 1 19851 : iKonigl & KMti3 
[T994 : Elvis 2000; Proga 2000, 2007). However, these are probably 
tenuous, with an initial mass-loading < Mbh (although in at least 
some cases, these outflows are extremely de nse, and might have 
much higher mass-lo ading factors; see e.g. iHall & HutsemekersI 
l2003l ; lHall et alj|2007b . It is not clear whether such "hot" outflows 
could efficiently entrain gas at larger radii. If most of the gas mass 
of the galaxy is at some appreciable fraction of the galaxy effective 
radius Re and in the form of cold, dense giant molecular clouds 
(GMCs), then it is difficult to imagine such a diffuse wind di- 
rectly "launching" the clouds out of the galaxy. It remains unclear 
whether, in fact, the momentum associated with the winds that are 
known to emanate from the central engi ne of a quasar is sufficient 
to unbind the cold ga s in th e host (see e.g.^Bau m & McCarthvl2OO0l; 



de Kool et al] 1200 ll [200 2; Steenbrugge et a ll |2005|; 



Holt et al.l 



2006; Gabel et al. 2006; Krongold et al. 2007a b; Batcheldor et al. 
2007 ; Tremonti et al. 2007; McKernan et al. 20()7; Gangulv et al., 
20071 ; IProchaska & Hennawill2009l) . 



In this paper, we argue that it is not necessary that the small- 
scale, high-velocity AGN outflows directly entrain any cold gas at 
scales ~ Re. Rather, so long as these are sufficient to drive a signif- 
icant wind in the "hot" diffuse ISM, then clouds will be effectively 
destroyed or deformed and "secondary" feedback mechanisms - 
namely the radiative effects of dust absorption and ionization - will 
be able to act efficiently on the cold gas at large scales. This will ef- 
fectively terminate star formation on a short timescale, with greatly 
reduced energy/momentum requirements for the "initial" outflow 
drivers. 



2 RADIATIVE FEEDBACK IN THE PRESENCE OF HOT 
OUTFLOWS 

Consider a typical galaxy, where the ISM gas is composed of a 
mix of diffuse warm/hot gas and cold clouds^ Radiation will al- 
ways act on the cold clouds (in the form of ionization and momen- 
tum injection from absorbed photons), but they may be too dense 
and self-shielding to be significantly affected. In such a case, one 
could invoke a blastwave or cold shell, driven by AGN feedback on 



^ By "diffuse" ISM, we refer to tfie warm/hot pliases of the ISM, witfi effec- 
tive temperatures T ~ lO'* — lO^K and densities nu ~ 10~^ — 10''cm~^. 
The "clouds," on the other hand, represent the cold phase (but not col- 
lapsing molecular cloud cores) with T < lOOK molecular gas and iih ~ 
10^ - lO'* cm . The volume filling factor of clouds is small, but by mass, 
they dominate the total gas mass; the hot gas with filling factor near unity 
has a typical mass fraction /hoi J; 0. 1 (see e.g. ,McKee & Ostrikeinig??.) . 



small scales, to entrain this material. Various models and simula- 
tions have shown that if feedback needs to directly launch a blast- 
wave in both the hot and cold gas together (sufficient to entrain 
most of the galactic gas), then an efficiency tie ~ 0.05 is the rel- 
evant value (for energy injection where the feedback E — rjEL = 
riEtrMBHC'; ~ 0.1 is the radiative efficiency). If the outflow is 
instead momentum-driven, the relevant value for driving hot-l-cold 
phases is 77p ~ 1 (p = iq^L/c). 

If, however, the "initial" feedback from the central source need 
only drive a wind in the low-density hot gas, and does not neces- 
sarily directly entrain the cold clouds, then the energy required will 
be much less. In both cases above, the efficiency needed simply 
scales linearly with the mass of material to be driven (scaling with 
its binding energy or momentum, respectively). So, if a fraction 
/hot of the gas is in the hot, diffuse ISM, and only that needs to be 
initially driven, we obtain 



77e = ^ ^ 0.05 /hot 



for energy-driven and 



L/c 



'/he 



(1) 



(2) 



for m omentum-driven outflows. For typical /hot ^0.1 (e.g. iBlitzl 
Il993h . this implies an order-of-magnitude reduction in the neces- 
sary feedback input for "interesting" behavior. And in detail, since 
the hot-phase gas is already virialized (rather than in e.g. a cold 
disk), the efficiency gains may be even higher. The question is then, 
if most of the mass is in cold clouds and only the hot gas is affected 
by this initial outflow, will any interesting behavior result? 

Consider a cold molecular cloud in a galaxy. The cloud has 
a mass Mr, and an initial ("equilibrium") effective spherical radius 
Ro- In the observed typical ISM, these are related by 

Me 300 Mq Rip, (3) 

(where Ro.pc = i^n / 1 pc; iBlitz & Rosolowskvll2006h . Since we will 
later allow the cloud to stretch or deform, define the instantaneous 
radius of the cloud as Re. The cloud resides at a spherical distance 
r from the center of a Hernquist ( 1990) profile bulge of total mass 
Mbiii and characteristic scale length Rcsi, which hosts a BH on the 
characteristic BH-h ost relations, Mbh ~ MEnMbui (^ibh ~ 0.0014; 
iHaring & Rixll2004h . The BH radiates with a luminosity 



^Qso = rhLEdd 



(4) 



where m is the dimensionless Eddington ratio and LEdd ~ 
1.3Mbh.8 X lO^'^ergs"' is the Eddington luminosity (Mbh.s ~ 
Mbh/10' Mq). 

We are interested in the case where the cold clouds are embed- 
ded in some kind of hot outflow generated by the "primary" AGN 
feedback. Specifically, assume that the quasar somehow succeeds 
in driving an outflow through the diffuse warm/hot ISM; to be con- 
servative, the outflow can be tenuous and we will assume that the 
outflow "impacting" the cloud carries negligible momentum com- 
pared to the binding momentum of the cloud. In other words, as- 
sume that at some smaller scale, a wind or outflow is generated 
sufficient to sweep up the tenuous, diffuse ISM at large radii, but 
insufficient to affect cold dense clouds, which contain most of the 
ISM gas mass. 

This problem of the survival of cold clouds in a post-shock hot 
medium is well-studied in the context of s tar formation and super- 
nova feedback (see e.g. iKlein et al]| 19941 and references therein). 
In general, the collision of a shock or wind with velocity v, with a 



© 0000 RAS, MNRAS 000, 000-000 



Multistage AGN Feedback 3 



Incident 
Quasar 
Radiation 




Low-density Gas 



Diffuse Outflow 



• t~2-5t3 ^ 



Figure 1. An illustration from a simple simulation of the effects of a hot outflow on radiation feedback. Left: Initial dense cloud in pressure equilibrium with 
the diffuse ISM. Ionization and momentum flux from the quasar is negligible, because the effective surface area for absorption is small. Center: At f = f,, 
an outflow generated in the hot/diffuse ISM hits the cloud. Ambient pressure drops in the wake (low-pressure regions form on the perpendicular sides of the 
cloud). A mix of instabilities shortly cause the cloud to mix in the perpendicular direction. Right: After a few cloud crossing times, the mixing has increased 
the effective cross section by a factor [Rc/Rq)- > 10 in the perpendicular direction. The cloud may now (with lower density and higher cross section to absorb 
quasar photon momentum) be vulnerable to radiative feedback, and will be accelerated to v ~ Vesc . 



cloud of initial characteristic (quasi-spherical) radius Ro and den- 
sity contrast x (ratio of cloud density to external medium den- 
sity X = "c/no) will launch secondary shocks within the cloud 
with velocity Vs/x^^^- This defines a "cloud crushing" timescale 
fee = x^^^ Ro/vs- In the simple case of a pure hydrodynamic strong 
shock, if tec is much less than the characteristic timescales for the 
density to change behind the shock and Vsjx^^' is comparable to 
or larger than the characteristic internal velocities of the cloud, 
then the cloud will be stretched and "shredded" by a combina- 
tion of Rayleigh-Taylor and Kelvin-Helm holtz instabilitie s on a 
timescale ^a few fa (see e.g.lKlein et alll 994: Xu & St on^l 19951 : 
iFragile et alj |2004| : lOrlando et al.l l2005l : lNakamura et idri2006h . 
Given the definitions above, for v., ~ Vese, fee is much less than the 
dynamical timescales of interest at all the spatial scales of interest 
(fee «C 10^ yr for all r > x^^^Ro - i-e. for clouds not in the very 
nuclear regions). 

The material being mixed off of the surface of the cloud from 
these instabilities expands into and mixes with the low-pressure 
zones created by the passage of the shock on the sides of the cloud. 
This leads to an effective net expansion of the cloud by a factor 
~ x'^^ in radius in the perpendicular shock direction. Eventually, 
despite the initial compression, reflection shocks lead to a expan- 
sion by a factor --^ 2 in the parallel shock direction, bringing the 
original cloud material into an effective density and pressure equi- 
librium with the external medium. The surface area of the cloud 
can increase dramatically; for our purposes, we are interested in 
the effective cross section the cloud presents to the perpendicular 
shock direction. The radii defined above should be thought of in 
this manner: the initial cloud cross section to the hot shock isttRI; 
post-shock, the effective cross section owing to this expansion and 
equilibration isivR^ ^ X'^Ro- 

We illustrate this behavior with a simple toy model system 
in Figure [T] Specifically, we show an example of a hydrody- 
namic simulation of an i dealized system, using the ZEUS code 
jStone & NormanI Il992al lbl). The initial conditions consist of a 
Plummet sphere cloud embedded in a uniform background, with 
density contrast of x = 100 (peak density of the cloud relative to 
background), with the initial system in pressure equilibrium (uni- 
form pressure), and periodic boundary conditions in a large grid. 
At time t = ts the low-density material is rapidly accelerated into 
a mach ~ 2 wind. Color encodes the gas density, from black (the 
arbitrary background density) to red (the initial maximum). Note 



that the example shown is purely for illustrative purposes - we do 
not include many possible complexities, such as gas cooling, star 
formation, or magnetic fields. The behavior of clouds in response 
to outflows with such sophistications has been extensively studied 
in the references above, and more detailed extensions of such sim- 
ulations to the regime of interest here will be the subject of future 
work. Nevertheless, this simple experiment illustrates much of the 
important qualitative behavior. 

The qualitative behavior we care about - the mix- 
ing/stretching/deformation of the cloud in the perpendicular shock 
direction leading to an increase in the effective cross section of the 
cloud - is in fact quite general. Simulations have shown that the 
same instabilities operate regardless of whether the "hot outflow" is 
a strong shock, weak shock, or wind (since we assume the hot ma- 
terial is being unbound in this wind, it cannot be substantially sub- 
sonic). The timescale of cloud expansion increases by a factor of a 
few in the weaker wind case, but it i s still much less than the rele- 
vant l ocal galactic dynamical times taein et al.ll 19941 : 1 Jones et al.1 
1 19961) . The process is also similar in the case of a cloud being 
impacted by AGN jets, despite the different densities, tempera- 
ture s, and magnetic field states associated with jets and "bubbles " 
(see lKrause & Alexanderll2007l : lAntonuccio-Delogu & SilklboOSh . 
A cloud could in principle be stabilized against suc h instabilities 
by being strongly magnetically domi nated (Mac Low et af]| 19941 : 
lJun & Jonedl 19991 : iFragile et al.ll2005l : IShin et al. 2008). However, 
in this limit, as the hot outflow sweeps up material, the pressure of 
the diffuse ISM trailing the outflow will decline as a steep func- 
tion of time ?/fcc dOstriker & McKeel 19881) . Since, in this limit, the 
cloud is then over-pressurized, it will expand isothermally as the 
exterior post-shock pressure drops (the free expansion/equilibration 
time of the cloud being short compared to the other timescales of 
interest). Because this stops when the system is equilibrated, the 
"effective" net expansion of Rc is the same as in the hydrodynamic 
shock case, even though the details are quite different. 

If nothing more were to happen to the cloud, this would only 
suppress star formation for a short time. The cooling instabilities 
that produced the cloud in the first place would operate. In the "typ- 
ical" ISM, clouds mix in the wake of stellar or supernovae-driven 
outflows until they reach equilibrium with the ISM and re-cool into 
new clouds. 

However, we are interested in all of this occuring in the 
background of a luminous AGN, which will both ionize and ex- 
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ert a radiation pressure force. Tiie cloud - especially a realistic 
cloud with a large dust mass and corresponding opacity - is op- 
tically thicl0 to the quasar radiation, with an effective cross-section 
A£7 ~ [ttR^) I (47rr^). There is therefore an inescapable deposition 
of photon momentum from the radiation field with a deposition 
rate piad = L-^hs/c where Labs = Lqso Af2. Comparing this to the 
gravitational force Fgrav = —Mcd(j>/dr defines an effective Edding- 
ton limit for the cloud: if the absorbed flux exceeds some limit the 
cloud will be unbound (equivalently, the absorbed momentum in a 
single dynamical time, over which the cloud could redistribute that 
momentum, will exceed the cloud binding momentum ~ M^Vesc). 
This limit is when the two are equal: equivalently 



Lqso ttRc 



■-Mc 



(5) 



Assuming that the cloud lies on the observed size-mass rela- 
tion (Equation [3} and that the galaxy lies on the Mbh — Mtui rela- 
tion, this reduces to the criterion for unbinding the cloud: 



RqJ ~ 



r + R,if 



(6) 



In other words, a cloud on the "normal" size mass relation at large 
radii r ~ i?eff is sufficiently dense and sufficiently high column- 
density to avoid being unbound by radiation pressure. But if the ef- 
fective size of the cloud (the effective coupling surface area) could 
be increased by a factor of a couple, or the effective column low- 
ered, the cloud would rapidly be unbound by the incident radiation 
field momentum. This condition is easily satisfied in post-shock 
clouds. 

Figure [T] also includes and illustrates this effect. Specifically, 
we include a very simple, time-independent momentum deposition 
rate in the "facing" cells to the incident radiation field, which we 
simply approximate as all cells with a density above > 3 times 
the initial background density but with no cell at x < x, above this 
density (i.e. implicitly assuming that such a cell would shield the 
cells "behind" it with respect to the incident radiation). The mag- 
nitude of this is initialized such that at / < ts the "total" deposition 
rate over the surface of the cloud is equal to a small fraction of the 
"binding" momentum over the dynamical time (assuming I'j ^ Vc), 
~ O.Ol McVs/{Rc/vs). But the details make little qualitative differ- 
ence to the global acceleration. 

A similar effect pertains to the ionization of the cloud (al- 
though this is not explicitly included in Figure[T](. Ignoring geomet- 
ric effects of photon diffusion, the volume of a cloud of mean den- 
sity ric ionized is Vi„„ —N/n^fi, where /3 ~ 2 x 10~ cm^ s~ ' is the 
recombination coefficient for gas at the temperature for hydrogen 
ionization (Te ~ 10"* K) and N is the rate at which ionizing photons 
hit the cloud. The total rate of production of ionizing photons from 
the quasar is Nq = XL/hvgn (A ~ 0.07 co mes from a proper inte - 
gration over the quasar spectrum; here from lHopkins et al 
and a fraction Afi are incident on the cloud. Together with the typ- 
ical values above, this implies that clouds will be ionized to a depth 



' For any incident spectrum where a significant fraction of the inci- 
dent energy is in the optical/UV or higher wavelengt hs, the effective op- 
tical depth from dust within the cloud will be r ~ 1 iMunav et ai]|2005l : 
[Thompson et al.ll2005h - this is simply a statement that the clouds will be 
optically thick to some portion of that SED (whether ionizing photons in 
the far UV, or, if the cloud is not yet ionized, then there is dust which will 
absorb in the optical and IR) and re-radiate that energy. The cloud can ef- 
fectively be thought of as a single absorbing "mega-grain" with effective 
cross section ~ -kRI. 



Rc 



Rc\^ 



RoJ 



(7) 



Give the cloud size-mass relation, this is equivalent to the statement 
that all clouds below a mass Mc < 10** Mq (Rc/Ro)"^" will be self- 
shielded at r > /?cff. For typical clouds, the depth ionized is clearly 
quite small; but there is a steep dependence on cloud radius. As Rc 
increases, the ionized depth increases by a factor oc R^ owing to the 
increased photon capture cross-section and a factor oc rI owing to 
the decreased density lowering the recombination rate. 



3 IMPLICATIONS IN A GLOBAL FEEDBACK 
SCENARIO 

There are many caveats to the simplified derivations above: clouds 
have some size and mass spectrum, and are distributed at various 
radii, with the background quasar changing in time. Nevertheless, 
embedding this in more detailed models for AGN feedback, the 
results are interesting. 

Consider an ~ L « bulge w i th Mb u\ = 10"Mq and Mbh = 
1.4 X 10** Mq, with a lHernquistl j 19901) density profile and scale 
radius R^k = 4kpc. Assume gas traces stars (initially), with mass 
fraction /gas = 0. 1, and that 90% of the gas is in cold clouds while 
10% is in a hot diffuse phase (which we assume is in hydrostatic 
equilibrium). This yields a density profile of hot or cold gas of 



P, =/,/gas 



Mbal 



Re 



(8) 



27r R{R + Rcy 

where represents the fraction in either the hot or the cold phase, 
(/ho. = 0.1, /cold = 0.9). 

At a given radius (within some small radial annulus), the 
cold gas (with mean volume density given above) is assumed to 
be locked into cold clouds with a small volume filling factor. 
The clouds are initially placed on the observed size-mass relation 
(Equation[3] determining an initial radius Ro for each cloud of mass 
Mc), and are distributed in mass according to the observed mass 
spectrum dN/dMc oc M7'^ (up to a maximum Mc = lO' Mq; see 
JRosolowsky 2005! ^rid references therein). Given some total cold 
gas mass pcoid per unit volume in an annulus, this mass spectrum, 
integrated from arbitrarily small cloud mass (it makes no differ- 
ence if we adopt some lower-mass cutoff) to the maximum Mc, 
must integrate to pcoid per unit density; i.e. this determines the num- 
ber density «c of clouds at each galaxy radius R and mass interval 
Mc Mc + dMc (and corresponding initial cloud radius Ro{Mc)). 

At a time ? = 0, we assume that the BH "turns on," radiating 
(initially) at the Eddington limit. We allow it to drive a shock/wind 
th rough the diffuse ISM a ccord ing to the analytic solutions derived 
in iHopkins & Hemguistl ( l2006h . In these models, the AGN is as- 
sumed to drive a simple, Sedov-Tay lor- type outflow via any "small- 
scale" feedback channel. The detailed behavior is derived and 
compa red with hydrodynamic simulations in Hopkins & Hernquis^ 
( l2006h . but can be simply summarized as follows: the BH, once 
on, couples a fraction 77 ~ 0.05 /hot of its luminous energy to the 
diffuse-phase gas in its vicinity. Because the spatial and timescales 
in the vicinity of the BH are small compared to the rest of the 
galaxy, this appears to the galaxy as a point-like energy injection in 
a hot medium. The result is therefore a roughly self-similar (power- 
law) Sedov-Taylor-type outflow. 

Figure [2] illustrates some of the basic behaviors of this out- 
flow. The energy injection leads to a shock that expands out- 
wards with radius Rs{t), in an approximately power-law fashion 
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as Rs oc t", where a is a function of the local density profile, gas 
equation of state, and is coupled (weakly) to the declining en- 
ergy injection rate of the BH (Figure |2l bottom left). For typical 
galax y density profiles, a 7 = 5/3 gas, and the conditions assumed 
here, iHopkins & Hemguistl ilOOdtl show Q « 4/5. In the wake of 
the expanding bubble/shock, the post-shock density drops in a re- 
lated power-law manner. Since the spherical accretion rate onto 
a BH scales roughly oc p (in Bondi-Hoyle accretion; or similarly 
oc E for viscous accretion from a disk), the accretion rate, and 
hence luminosity L = 0. 1 Mbh c^, will decay as well (Figure [2l top 
left). Again, we refer to the full derivation in lHopkins & Hemguistl 
I2OO6 ) for details, but the self-consistent solution derived therein 
can be approximated as L oc [1 + t/tg]'^^, with rjL ~ 1.5 and 
tQ ^ 1 — 5 X 10'' yr for the parame ters here (con s istent with var- 
ious observational constraints; see iMartinillBoi ; IYu et alJllOO^ ; 
iHopkins & Hemquistllioogh . The value of rjt follows from a and 
generic behavior of Sedov-Taylor post-shock gas, and tg is simply 
related, modulo appropriate numerical coefficients, to the local dy- 
namical time near the BH radius of influence. 

The shock therefore crosses a radius r at a time (when 
Rs — r). In the wake of the shock, the post-shock ambient pressure 
Pext will drop, reflecting the density decline from material being 
blown out (Figure|2l top center). Again, this approximately follows 
a standard decline in the wake of a Sedov-Taylor blastwave; the ex- 
act solution for the pressure internal to the blastwave unde r the con- 
dition s here must be obtained numerically, but Ostrike r & McKeel 
( Il988h show that it can be approximated as a double power-law. 
Roughly speaking, there is a rapid drop in pressure in the imme- 
diate post-shock region, as the thin shell at the front of the blast- 
wave clears the diffuse material away from the region, and the 
pressure declines as a steep power P oc {[t — ts\/tc)~^ , where tc 
is the crossing time of the shock relative to the cloud (~ Rc/vs) 
and /3 ~ 3 — 5 (the exact index depends on the local density profile 
slope and rate of decay of the driving source, so is not the same at 
all radii). This is followed by a more gradual decline, as the diffuse 
medium internal to the shock relaxes, is heated, and expands, with 
P oc ([f — ts\/ts)~^ and /3' ~ 2 (again, for the detailed derivation 
of the double power-law structure in the wake of the blastwave, for 
the conditi ons of the feedback-driven b lastwaves considered here, 
we refer to lHopkins & HemguistlllOO^ . 

In the wake of the shock, with a rapidly declining background 
pressure and density, the cloud will be mixed and effectively in- 
crease its surface area. We can solve for the behavior of each cloud 
- at least the key parameter of interest, the effective radius of the 
cloud in the direction perpendicular to the sh ock {Rc) as a func- 
tion of time, according to the approximations in lKIein et al.l ( 1 19941) 
in the wake of the shock defined above. If the cloud could some- 
how resist shredding (say, via sufficient magnetic field or turbulent 
support) this would be trivial: the cloud (initial pressure equilib- 
rium) cloud would expand isothermally (e.g. conserving total mag- 
netic energy) as it is now over-pressurized, such that pressure equi- 
librium would be conserved. For a background pressure declining 
with some power law P oc {t/ts)~^, this implies expansion of the 
cloud with Rc/Ro = [t/tsf^ (generally, Rc/Ro = {Pc/Po)~^^\ for 
isothermal expansion). If the cloud were supported by thermal en- 
ergy with no new inputs this index would be slightly modified (for 
adiabatic expansion, Rc/Ro ~ (-Pc/Po)'^^) but the behavior is qual- 
itatively similar. Technically this approximation assumes that the 
time for the cloud to equilibrate is short relative to the timescale on 
which the background is changing, but this is easily satisfied. The 
cloud expands/equilibrates on its internal crossing time, given by 
the effective sound speed as ~ Rc/cs.^fs; for quasi-virial clouds this 



is simply the dynamical time 1 l\JG p where the effective density p 
follows from the observed size-mass relation (Equation [3j- Using 
the observed values, this gives a timescale of ^ 0.1 — 1 x 10^ yr 
(for cloud sizes ~ 0.1 — lOpc). Compare this to the character- 
istic timescale for the evolution of the background, a few f, (it- 
self of order fdyn, the galaxy dynamical time at R). For a typical 
a ~ 200kms~' spheroid, ts 2> Rcjcs.ai for all radii > lOOpc - 
in other words, this condition is easily satisfied at the radii ~ Re, 
which contain most of the mass of the galaxy. Figure [2] (fco/tom 
center) shows how the clouds will expand in effective cross sec- 
tion (~ P^e)-, relative to their initial sizes, given this declining back- 
ground pressure for the simple isothermal case. 

For the more complex case of cloud shredding, it turns out 
that, in aggregate, a similar scaling obtains. The characteristic time 
for the cloud to effectively be mixed via instabilities and so effec- 
tively increase its cross section is a few cloud-crushing timescales 
tec = y^^'^Rfilvs. But since the shock velocity is of order the galaxy 
escape velocity for the "interesting" diffuse outflows considered 
here (v, ~ a few a), and typical a ~ 200 kms^', compared with a 
typical effective sound speed of a virialized cloud ~ 1 — lOkms"', 
this time is almost always much shorter than (or at least compa - 
rable to) the cloud dynamical time. Following iKlein et al.l ( Il994h . 
cloud shredding will equilibrate when the system expands by a fac- 
tor ~ x'^^ in the perpendicular direction (and a small, ~ constant 
factor in the parallel direction), where x the initial density con- 
trast; in other words, until the effective density and pressure drop to 
approximate equilibrium with the background. Thus, for timescales 
~ f, over which the background is evolving, long compared to the 
cloud-crushing time, we can consider the systems to effectively ex- 
pand with an average effective radius scaling in the same way in 
equilibrium with the external/hot medium background pressure (i.e. 
similar effective net expansion, averaged over these timescales, as 
in the isothermal expansion case). 

We then solve for the behavior of each cloud in the wake of 
this h ot outflow, according to the approximations in iKlein et al.l 
\mA) and §[2] (Figure H top right). In particular, given the time- 
evolution in Rc/Ro shown in Figure |2l we use the scalings derived 
in § [2] to estimate the fraction of the cloud (at some initial radius 
r) which will be ionized (i.e. fan ~ 'lionized /^c from Equation |7J 
where the AGN accretion rate m and cloud expansion Rc / Ro are 
given as a function of time above, the initial cloud radius r is one 
of those specified in the Figure, and we chose a representative ini- 
tial cloud with radius Ro.pc ~ 1 for illustrative purposes). We also 
show the relative strength of radiation pressure on the cloud, i.e. the 
radiation pressure relative to the local Eddington limit (that which 
would unbind the cloud), P,ad/Edd = 0. 14ra {Rc/Roy l{r + Re) /r]' 
(re-arranging Equations|5]&[6j where again rn, Rc/Ro, and r for the 
clouds is given). Since these both scale steeply with the increas- 
ing effective cloud size (oc R^ and R^, respectively), both increase 
rapidly in time. At early times, only the clouds within a narrow re- 
gion ~ lOOpc around the quasar are efficiently ionized, and the ef- 
fects of radiation pressure are we ak - similar to what is observed in 
the n arrow-line regions of AGN jCrenshaw et"ai]|2000l ; | Rice et al.l 
I2OO6I) . This changes rapidly in the wake of the outflow at large radii 
- the deformation induced makes the clouds vulnerable to ioniza- 
tion and radiative momentum driving. 

Integrating over the entire cloud population and galaxy mass, 
i.e. integrating over the initial cloud mass or size (Ro) spectrum at 
each galactic radius r, and then over the total cold gas density at 
each radius r, we obtain the total fraction of cloud mass that can be 
ionized or effectively accelerated by radiation pressure (we define 
the latter as the integral of mass in clouds where Frad/Edd > 1). 
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Figure 2. Effects of a two-stage feedback model in the wake of a quasar-induced hot outflow, with parameters of a typical ~ L, quasar. Top Left: Quasar 
Eddington ratio as a function of time, after feedback from the quasar begins to drive an outflow in the diffuse ISM (time t = 0). Bottom Left: Radius of the 
hot/diffuse outflow as a function of time (compare the galaxy effective radius i?eff). Top Center: Change in pressure of the (post-shock) diffuse ISM at a given 
radius r from the BH, at time t /is (Is is the time when the hot outflow first reaches r). We show three different radii: r = lOOpc, 1 kpc, and lOkpc. Bottom 
Center: Change in effective cross-section of a typical {Rq ~ Ipc) cloud at each r versus time t/ts- Top Right: Ratio of the radiation pressure force to the 
Eddington limit for expulsion of each cloud (thick; P„(]/Edd), and fraction of the cloud ionized (thin; /io„). Bottom Right: Total fraction of cloud mass that 
can be expelled by radiation pressure or ionized (we consider each separately), as a function of time. The curves reflect the integral over the observed mass 
spectrum of clouds at each radius, integrated over all radii following the galaxy density profile. 



Figure [2] {bottom right) shows this as a function of time given the 
above outflow properties. Most of the cold mass can be effectively 
accelerated by radiation pressure at large radii in a timescale ~a 
few galaxy dynamical times, despite the declining AGN luminosity. 
We have experimented with varying the exact parameters adopted 
here, and find that the qualitative results are robust. In a couple of 
dynamical times, ~ 90% of the original cold cloud mass becomes 
vulnerable to secondary radiative feedback; i.e. these numbers - the 
fraction that can be ionized and/or the strength of radiation pressure 
approach or exceed unity. 

This and the previous results also provide an important check 
of an implicit assumption in this model: that the cloud acceler- 
ation time is long relative to the timescale for the clouds to ex- 
pand/equilibrate. If this were not true, the two behaviors could not 
be treated independently, and the physical consequences are un- 
clear (it is possible, for example, that each "parcel" of the cloud 
which is mixed or stripped off by instabilities would rapidly be ac- 
celerated, leading to the edges of the clouds being effectively blown 
away or stripped but giving little acceleration to the cloud core). As 
noted above, the expansion/equilibration time is given by the cloud- 
crushing time, comparable to the internal cloud crossing/dynamical 
times < 10* yr. In comparison, the acceleration times are of order 
a couple to a few the dynamical time at r in the galaxy, which 
(as we discuss above) are generically much larger than the cloud 
crossing time at all radii ~ Re (where most of the galaxy mass is 
located), indeed all radii > lOOpc (i.e. all radii which are not al- 
ready affected by feedback even without a diffuse outflow). In a 
global sense, most of the mass is accelerated and/or ionized over 
a timescale 10^ — 10^ yr, much longer than the crossing/collapse 
times of all but the most massive molecular cloud complexes. 

In fact, this acceleration time is a relatively long time, at large 



scales (~a few 10** yr), and the AGN luminosity has correspond- 
ingly decayed to ~ 1% of the Eddington limit. The model above 
accounts for this, but an important question remains if real AGN 
can sustain even this level of energetic output over these time inter- 
vals. If, for example, AGN switch to a radiatively inefficient state 
above or around this accretion rate, the driving will suddenly van- 
ish. Clearly, this is an interesting regime; better knowledge of how 
feedback-induced hot outflows and subsequent lightcurve evolution 
proceed will be important to understanding both how dramatic the 
effect on the galaxy will be and, potentially, how much variation 
there may be between galaxies. 



4 DISCUSSION 

"Feedback" from bright AGN is a topic of fundamental interest for 
galaxy evolution, but it remains unknown whether or not any of the 
obvious candidate feedback mechanisms are capable of effectively 
coupling to cold molecular gas, especially at kpc scales, the domi- 
nant reservoir for star formation. Here, we demonstrate that it is at 
least possible that the cold gas reservoir is destroyed and/or blown 
out of the galaxy despite inefficient coupling of "initial" feedback 
mechanisms that originate near the BH. 

If some coupling of energy or momentum near the BH 
- whether from e.g. Compton heating, radiation pressure, BAL 
winds, jets, or resonant line-driving - can generate a wind or 
shock/blastwave in the warm/hot ISM, then when the outflow 
passes by a cold cloud, even if it does not directly entrain the ma- 
terial, it will generate various instabilities that "shred" the cloud 
and mix it, efficiently enhancing the cloud cross section in the per- 
pendicular direction. Even if the cloud is magnetically supported or 
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extremely dense and able to resist instabilities, there is still a grow- 
ing pressure imbalance that drives the cloud to expand in the same 
manner. 

This is well-studied in the context of supernovae-driven 
winds, but there is an important difference in the presence of a 
bright quasar. The effective increase in cross section means that 
momentum driving and ionization heating from the quasar radia- 
tion is quickly able to act in much more dramatic fashion on clouds 
that were once too dense and too small (or at too large a distance 
from the black hole) to be perturbed by the radiation field. This 
effect can have dramatic implications for star formation in quasar 
host galaxies. 

Because radiation pressure always acts, this means the energy 
needed in "initial" feedback from the central source to e.g. drive 
winds in the low-density hot gas will be much less than if it were 
expected to act directly on the cold clouds. We show that the en- 
ergetic or momentum driving requirements for the initially driven 
feedback are reduced by at least factor /hot ~ 0. 1 (the mass frac- 
tion in the hot diffuse ISM); i.e. rather than the canonical ~ 5% of 
the radiant energy (~ 100% momentum) needed in the initial out- 
flow if it were to entrain the entire gas supply directly, only ~ 0.5% 
(~ 10% momentum) is sufficient to drive the hot gas and enter the 
regime of interest here. Another way of stating this is, for accretion 
with an Eddington ratio m and BH mass Mbh relative to the expec- 
tation (Mbh) from the Mbh — o relation, the relevant outflows will 
be driven (and star formation suppressed) when 

where 77 is the feedback efficiency {E = r/L). Given this criterion, 
that there is sufficient momentum in photons for the "secondary" 
feedback to act is gauranteed for all but the most extremely gas- 
dominated systems. 

Note that the derivation here pertains to large clouds (Ro > pc), 
observed to be in rough pressure equilibrium with the ambi- 
ent medium and containing most of the ISM mass. Dense cores 
(Ro <^ pc) are observed to be in self-gravitating collapse; these will 
continue to collapse and form stars on a very short timescale de- 
spite a diffuse outflow. The important thing is that no new cold gas 
reservoir of large clouds will be available to form new cores. 

We have also neglected the possibility that galaxies are highly 
self-shielding. For example, in dense nuclear star-forming regions 
in e.g. ULIRGs, the column densities are so high (A'h > 10^^ cm~^, 
see e.g. Komossa et al.. 2003 ; Li et al..2007) that the quasar can do 
little until star formation exhausts more of the gas supply. In disk- 
dominated galaxies, gas near R^f can be similarly self-shielded. If 
the radiation is isotropic, then for some disk mass and gas frac- 
tion, only a fraction ~ h/R (the fractional scale height at R) will 
couple to the relevant area (and the radiation may, in fact, be pref- 
erentially polar, yielding even lower efficiency). Only the most gas- 
poor disks, or the central regions of disks (where systems are typi- 
cally bulge-dominated) will be affected by the coupling efficiencies 
above. 

What we outline here is a simple model for the qualitative 
physical effects that may happen when cold clouds in the ISM en- 
counter a hot outflow driven by an AGN. More detailed conclu- 
sions will require study in hydrodynamic simulations which incor- 
porate gas phase structure, cooling, turbulence, self-gravity, radi- 
ation transport, and possibly (if they provide significant pressure 
support) magnetic fields. Detailed effects which we cannot follow 
analytically, such as e.g. self-shielding within thin, dense fingers in 
Rayleigh-Taylor or Kelvin-Helmholtz instabilities may alter the ef- 



fects of the radiation field on the cloud materia and change our con- 
clusions. Nevertheless, our simple calculations here demonstrate 
that the process of cloud deformation in the wake of a hot out- 
flow can have dramatic implications for the susceptibility of those 
clouds to other modes of feedback, and should motivate further 
study. 

According to these simple considerations, outflows driven by 
AGN feedback may in fact be "multi-stage" or "two-tiered", with 
an initial hot Shockwave or strong wind driven by feedback mech- 
anisms near the BH, which is then supplemented by a succes- 
sive wind driven out as clouds in the wake of the former are de- 
formed/mixed and increase their effective cross-section to the AGN 
luminosity. The characteristic velocity of this secondary outflow, 
which will carry most of the mass, should be ~ Vcsc at the radii 
of launching (~ lO^kms^'), and it will behave similarly to out- 
flows from star formation. Indeed, because the driving occurs at 
large radii, it is not clear whether it could be distinguished from 
stellar-driven outflows at all, except indirectly (e.g. in cases where 
the observed star formation is insufficient to power the outflow). 
Charactersitic timescales are ^ a few /dyn of the galaxy, so much 
of the outflow occurs as sub-Eddington luminosities (as the AGN 
fades in the wake of launching the "primary" outflow) and the sys- 
tems will not appear gas-depleted until they have evolved by a sig- 
nificant amount (~few 10*yr, at Eddington ratios ~ 0.01 typical of 
"quiescent" ellipticals). These processes should nevertheless imply 
effective shutdown of star formation and destruction/heating of the 
cold gas supply in "massive" BH systems - bulge-dominated sys- 
tems on the Mbh — o" relation that have recently been excited to 
near-Eddington luminosities. 

Most intriguing, this reduces the energetic requirements for 
the "initial" feedback - whatever might drive an outflow in the 
hot gas from the vicinity of the BH - by an order of magnitude. 
Our estimates suggest that coupling only a fraction ~ 10^^^ of 
the luminosity of the AGN on small scales would be sufficient to 
drive such a hot outflow and then allow ~ 100% of the radiative 
energy/momentum to couple to cold gas. If, for example, quasar 
accretion-disk (or broad-line) winds (with characteristic velocities 
V ~ 10"'kms~') do not immediately dissipate all their energy, then 
the hot outflows we invoke would be generated with a mass-loading 
in such winds of just ~ 0. 1 Mbh, a fraction of the BH accretion rate. 
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